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We  question  recent  studies  invoking  the  existence  of  a  traditional  “logarithmic 
surface  layer”,  or  log  layer,  in  the  boundary  layer  of  the  rapidly-rotating  core 
of  a  hurricane.  One  such  study  argues  that  boundary-layer  parameterization 
schemes  that  do  not  include  a  log  layer  are  “badly  flawed”.  Another  study 
assumes  the  existence  of  a  log-layer  to  infer  drag  coefficients  at  hurricane  wind 
speeds.  We  provide  theoretical  reasoning  supported  by  observational  evidence 
as  to  why  significant  departures  from  the  normally-assumed  logarithmic  layer 
might  be  expected,  questioning  its  use  in  the  inference  of  drag  coefficient  at  high 
wind  speeds  and  laying  bare  suggestions  that  hurricane  models  using  boundary- 
layer  schemes  that  do  not  represent  the  log  layer  should  not  be  used.  The 
ramifications  of  these  findings  for  hurricane  modelling  are  discussed.  Finally, 
we  draw  attention  to  a  study  examining  a  range  of  boundary-layer  schemes 
demonstrating  that  a  recently  articulated  boundary-layer  spin-up  mechanism 
transcends  the  presence  of  a  log  layer. 
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1.  Introduction 

The  importance  of  the  boundary  layer  in  tropical  cyclones 
has  been  recognized  for  several  decades  because  the 
frictional  breakdown  of  gradient  wind  balance  leads  to 
strong  inflow  in  the  layer.  This  inflow  converges  moisture 
evaporated  from  the  sea  surface  to  feed  the  deep  convective 
clouds  in  the  storm’s  inner  core. 

The  boundary  layer  is  a  key  element  of  Ooyama’s 
seminal  axisymmetric  tropical-cyclone  model  (Ooyama 
1969),  Carrier’s  hurricane  model  (Carrier  et  al.  1971, 
Carrier  1971a, b)  as  well  as  Emanuel’s  (1986)  model  for 
an  axially-symmetric,  steady-state  hurricane.  Emanuel’s 
model  became  the  basis  for  a  widely-used  theory  for  the 
potential  intensity  (PI)  of  a  hurricane,  i.e.  the  maximum 
gradient  wind  speed  (Bister  and  Emanuel  1998,  Emanuel 
and  Rotunno  2011).  In  all  of  these  models,  the  boundary 
layer  is  treated  as  a  layer  of  air  of  constant  density.  In 
both  Ooyama’s  and  Emanuel’s  models  the  boundary  layer 


is  treated  as  a  layer  of  constant  depth  and  with  vertically- 
uniform  properties  and  explicitly  (tacitly)  assumes  that  the 
layer  is  in  gradient  wind  balance  (Smith  et  al.  2008). 

Recently,  Smith  et  al.  (2009)  have  demonstrated  that 
the  role  of  the  boundary  layer  extends  beyond  that  of 
converging  moisture:  it  has  a  dynamical  role  in  converging 
absolute  angular  momentum1,  M.  Although  M  is  not 
materially  conserved  in  the  boundary  layer,  large  tangential 
wind  speeds  can  be  achieved  there  if  the  radial  inflow  is 
sufficiently  large  to  bring  the  air  parcels  to  small  radii  with 
minimal  loss  of  M.  This  spin  up  mechanism,  while  coupled 
to  the  interior  flow  via  the  radial  pressure  gradient  at  the  top 
of  the  boundary  layer,  is  tied  fundamentally  to  the  dynamics 
of  the  boundary  layer,  where  the  flow  is  not  in  gradient  wind 
balance  over  a  substantial  radial  span.  It  was  shown  that  this 


1M  is  defined  as  rv  +  I  fr 2 ,  where  r  is  the  radius,  v  is  the  (azimuthally- 

averaged,  storm-relative)  tangential  wind  speed  and  /  is  the  Coriolis 
parameter. 
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mechanism  accounts  for  the  occurrence  of  the  maximum 
tangential  wind  in  the  boundary  layer,  a  feature  that  has 
been  found  also  in  observational  studies  (Montgomery  et  al. 
2006,  Kepert  2006a,b,  Bell  and  Montgomery  2008,  Sanger 
et  al.  2012). 

The  idealized  numerical  calculations  of  Smith  et  al. 
(2009)  employed  a  relatively  simple  bulk  boundary-layer 
parameterization  scheme,  albeit  more  sophisticated  than 
Emanuel’s  slab  model.  For  this  reason,  the  calculations 
were  repeated  by  Smith  and  Thomsen  (2010)  using  a 
range  of  boundary  layer  schemes  having  various  degrees  of 
sophistication.  While  the  latter  study  showed  quantitative 
differences  in  the  intensification  rate,  mature  intensity  and 
certain  flow  features  in  the  boundary  layer  for  different 
schemes,  in  all  cases  the  maximum  tangential  wind  was 
found  to  occur  close  to  the  top  of  the  inflow  layer  (see  Smith 
and  Thomsen,  Figure  2  and  Figure  1  herein),  implying 
that  the  boundary-layer  spin-up  mechanism  articulated  by 
Smith  et  al.  (2009)  is  robust  and  not  dependent  on  a 
particular  scheme.  Similar  results  were  obtained  by  Braun 
and  Tao  (2000)  and  Nolan  et  al.  (2009a, b)  in  case  studies 
of  two  particular  hurricanes,  where  different  boundary-layer 
schemes  were  compared.  While  a  range  of  schemes  were 
investigated  in  all  these  studies,  some  relatively  crude  and 
others  rather  sophisticated,  none  of  the  studies  went  so  far 
as  recommending  a  particular  scheme. 

In  an  effort  to  address  this  issue,  Kepert  (2012)  compared 
a  range  of  boundary-layer  parameterization  schemes  in  the 
framework  of  a  steady-state,  height-resolving,  boundary- 
layer  model  in  which  the  tangential  wind  speed  at  the  top 
of  the  boundary  layer  is  prescribed  and  assumed  to  be  in 
gradient  wind  balance.  One  outcome  of  his  study  as  stated 
in  his  abstract  is  that  “  ...  one  popular  class  of  schemes 
is  shown  to  be  badly  flawed  in  that  it  incorrectly  predicts 
the  near-surface  wind  profile,  and  therefore  should  not  be 
used.  Another  is  shown  to  be  sensitive  to  diagnosis  of  the 
boundary-layer  depth,  a  difficult  problem  in  the  core  of 
the  tropical  cyclone,  and  caution  is  advised.  The  Fouis 
boundary-layer  scheme  and  a  higher-order  closure  scheme 
are,  so  far  as  we  can  discern,  without  major  problems, 
and  are  recommended.”  In  his  conclusions,  Kepert  states 
that  “one  class  of  schemes,  representing  the  Bulk  and  Hi- 
Res  parameterizations2  available  within  MM53,  produces 
the  strongest  surface  inflow,  strongest  supergradient  jet, 
and  fails  to  produce  the  observed  near-surface  logarithmic 
layer”  and  “these  features  are  due  to  the  diffusivity  being 
a  maximum  at  the  lowest  model  level,  which  in  turn 
is  due  to  an  incorrect  parameterization  of  the  mixing 
length.  These  schemes  are  therefore  significantly  in  error  on 
observational  and  theoretical  grounds”  and  that  “  ...  it  would 
seem  prudent  that  such  studies  be  repeated  with  a  more 
reasonable  parameterization.”  Kepert  does  not  elaborate  on 
what  constitutes  “significantly  in  error”  (presumably,  this 
remark  applies  also  to  Emanuel’s  widely-used  PI  theory  that 
assumes  a  slab  boundary  layer)  and  we  question  here  the 
“observational  and  theoretical  grounds”  that  underpin  his 
claim. 

The  main  basis  of  Kepert’s  critique  of  many  schemes  is 
that  the  log-layer  has  to  be  satisfied  to  avoid  “significant 
error”  and,  for  consistency  with  a  constant  stress  layer, 
the  associated  mixing  length  must  increase  linearly  with 


2See  Kepert’s  paper  for  a  more  detailed  description  of  these  schemes. 

3MM5  refers  to  the  Pennsylvania  State  University/National  Center  for 
Atmospheric  Research  mesoscale  model. 
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depth  near  the  surface.  This  assertion  appears  to  be  founded 
on  an  observational  study  of  the  hurricane  boundary  layer 
by  Powell  et  al.  (2003)  and  on  laboratory  measurements 
in  non-rotating  boundary  layers  in  a  turbulent  fluid  (Von- 
Karman  1921,  Schlichting,  1979,  see  also  Stull  1988, 
Garratt  1992).  Using  a  composite  analysis  of  a  large  number 
of  Global  Positioning  System  (GPS)  dropwindsonde 
soundings  in  the  inner  core  of  storms,  Powell  showed  that 
the  logarithmic  layer  provides  an  acceptable  fit  to  the  wind 
speed  data  below  about  200  m  (see  his  Figure  1),  although 
there  is  a  large  scatter  in  the  wind  speed  data  and  Powell 
showed  only  data  points  at  each  height  and  not  individual 
vertical  profiles.  The  existence  of  such  a  layer  is  used  by 
both  Powell  (2003)  and  Holthuijsen  et  al.  (2012)  as  a  basis 
for  estimating  the  drag  coefficient  at  major  hurricane  wind 
speeds. 

We  are  unconvinced  by  this  “observational  support”  for 
the  ubiquity  of  a  log  layer  in  the  core  region  of  a  rapidly- 
rotating  vortex  for  several  reasons  articulated  below.  We  are 
unconvinced  also  by  the  theoretical  support  for  the  log  layer 
in  tropical  cyclones  asserted  by  Powell  et  al.,  Kepert  and 
others,  which  is  based  on  dimensional  analysis  and  assumes 
horizontal  homogeniety.  The  purpose  of  this  article  is  to 
revisit  the  interpretations  of  Powell  et  al.  (2003)  regarding 
the  log-layer  and  to  question  some  of  the  scientific  issues 
raised  by  Kepert’s  (2012)  study. 

The  structure  of  the  paper  is  as  follows.  In  Section  2 
we  review  the  derivation  of  the  log  layer  and  explain  why 
it  may  be  inapplicable  in  a  rapidly-rotating  vortex  and  in 
section  3  we  show  examples  of  inner-core  dropwindsonde 
soundings  that  do  not  support  the  existence  of  a  log  layer. 
Section  4  considers  the  implications  of  the  issues  raised 
for  modelling  the  hurricane  boundary  layer  and  section  5 
presents  the  conclusions. 

2.  Theoretical  considerations 

2.1.  The  log-layer  revisited 

The  derivation  of  the  log  layer  for  the  atmospheric  boundary 
layer  is  reviewed  in  a  classical  paper  by  Tennekes  (1973) 
and  is  based  on  an  asymptotic  similarity  theory  expounded 
by  Blackadar  and  Tennekes  (1968).  The  starting  point 
is  the  equations  of  motion  for  a  stationary,  horizontally- 
homogeneous,  barotropic  boundary  layer  flow  with  constant 
density  p,  which  is  forced  by  a  geostrophic  flow,  ug: 

-  /(«  ~Vg)  =  -T-(-u'w'),  (1) 

dz 

f(u  -  Ug)  =  (2) 

where  u  and  v  are  the  standard  Reynolds’ -averaged  zonal 
and  meridional  wind  components  in  the  boundary  layer,  ug 
and  vg  are  the  corresponding  geostrophic  wind  components 
at  the  top  of  the  boundary  layer,  /  is  the  Coriolis  parameter, 
and  0  is  the  height  above  the  surface.  The  expressions 
— u'w '  and  —v'w1  are  the  vertical  turbulent  momentum 
fluxes  of  zonal  and  meridional  momentum,  respectively 
(primes  denote  a  departure  from  the  mean  flow,  w'  being 
the  perturbation  of  vertical  velocity).  Taking  the  magnitude 

of  the  geostrophic  wind  as  G  =  +  vg,  the  surface 

roughness  length  as  z0,  and  the  surface  stress  pul  (here  u* 
is  the  surface  friction  velocity),  it  is  possible  to  establish 
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a  relationship  between  the  two  nondimensional  quantities: 
u*/G  and  Ro  =  G/(fz0),  the  surface  Rossby  number. 
Typically  Ro  »  1. 

Tennekes  (1973)  notes  that  these  equations  admit  two 
kinds  of  self-similar  solutions: 

1  zf/u*  finite  (i.e.  finite  relative  height  in  the  boundary 
layer),  but  with  z/ za  — >  oo  and  Ro  — >  oo.  Then,  to  a 
first  approximation,  the  wind  profile  is  asymptotically 
independent  of  Ro  provided  it  is  plotted  as 


(3) 


where  Fx  and  Fy  are  some  universal  functions  to  be 
determined.  This  is  the  scaling  for  the  part  of  the 
boundary  layer  above  the  surface  layer. 

2  z/z0  finite,  but  with  zf/u*  — >  0  and  Ro  — >  oo. 
Again,  to  a  first  approximation,  the  wind  profile 
is  asymptotically  independent  of  Ro  provided  it  is 
plotted  as 


u* 


u* 


(4) 


where  Fs  is  another  universal  function  to  be 
determined  and  it  has  been  assumed  that  the  surface 
stress  points  in  the  .x'-dircction.  This  is  the  scaling  for 
the  surface  layer. 

Tennekes  notes  that,  although  (3)  is  valid  only  well  outside 
the  surface  layer  (z/z0  — >  oo)  and  (4)  is  valid  only  inside 
the  surface  layer  {z/z0  finite),  they  must  have  a  region 
of  common  validity  if  Ro  is  large  enough.  This  region 
of  overlap,  in  which  zj za  — >  oo  and  zf/u*  —>  0  is  called 
the  matching  layer,  or  inertial  sublayer.  In  this  layer,  Eqs. 
(3)  and  (4)  and  all  their  derivatives  have  to  agree  with 
each  other.  Blackadar  and  Tennekes  (1968)  showed  that  the 
matching  is  possible  only  if  the  wind  profile  is  logarithmic 
with  height  and  that,  if  the  coordinate  axes  are  chosen  so 
that  the  surface  stress  is  in  the  x-direction,  Eqs.  (3)  and  (4) 
have  the  forms: 


U  —  Ug 


M* 


(5) 


and 


(6) 


where  B  and  k  are  constants,  the  latter  being  the  Von 
Karman  constant. 

In  the  case  of  a  steady  axisymmetric  vortex  with 
tangential  wind  speed  vg  (r)  in  gradient  wind  balance  above 
the  boundary  layer,  the  equations  analogous  to  (1)  and  (2) 
are 


_du  .  _  du 


_ v% 

Wdz  r 


qZ^U'W ')>  (7) 


u^t  +  ^  +  {f  +  Ca)u  =  Yz{-v'w^  (8) 

where  now  u  and  v  are  the  radial  and  tangential 
components  of  the  Reynolds’-averaged  wind,  respectively, 
and  v  has  been  replaced  by  the  agradient  wind  va  =  v 
vg.  The  derivation  is  assumes  the  normal  boundary-layer 


approximation  in  which  the  horozontal  pressure  gradient 
is  uniform  across  the  boundary  layer.  Now,  in  the  steady- 
state  case,  the  radial  flow  at  the  top  of  the  boundary  layer  is 
zero4  (i.e.  ug  =  0)  while  the  tangential  flow  is  a  function 
of  radius,  r.  Also,  £g  =  (1  /r)d(rvg)/dr  is  the  vertical 
component  of  relative  vorticity  of  the  gradient  wind. 

Although  a  scale  analysis  shows  that  the  nonlinear  terms 
in  Equations  (7)  and  (8)  cannot  be  neglected  in  the  boundary 
layer  of  a  tropical-cyclone  strength  vortex  (e.g.  Smith  1968, 
Carrier  1971a,  Vogl  and  Smith  2009),  let  us  suppose  for 
the  sake  of  the  current  discussion  that  they  can.  Then 
the  equations  are  similar  in  structure  to  (1)  and  (4),  but 
the  presence  of  the  radially-variable  coefficients  involving 
the  absolute  angular  velocity,  /  +  2 vg/r,  and  the  absolute 
vorticity,  /  +  £ g ,  of  the  gradient  wind  invalidates  the  scaling 
analysis  discussed  above.  In  a  rapidly-rotating  vortex,  these 
terms  are  dominated  by  the  contributions  2 vg/r  and  £g, 
which  typically  are  unequal  and  have  a  strong  radial 
variation,  except  possibly  close  to  the  centre,  where  the  flow 
may  be  in  approximate  solid  body  rotation.  In  this  case, 
the  assumption  of  horizontal  homogeneity  in  the  analysis 
of  the  equations  (1)  and  (2)  for  constant  /  is  no  longer 
valid.  Hence,  even  if  the  nonlinear  terms  are  ignored,  it  is 
by  no  means  obvious  to  us  that  a  similar  scaling  analysis 
can  be  applied,  since  additional  scales  including  the  radius 
and  relative  vorticity  of  the  gradient  flow  have  emerged. 

Alternative  derivations  of  the  logarithmic  velocity  profile 
in  the  surface  layer  are  common  in  the  literature  for  the 
case  of  a  homogeneous  flow  on  an  /-plane  (see  e.g. 
Panofsky  and  Dutton,  1984,  Stull  1988,  Garrett  1992, 
McWilliams  2006).  These  derivations  are  based  on  a  scale 
analyses  of  the  surface  layer,  itself,  without  considering 
a  formal  matching  to  the  boundary  layer  above,  although 
both  Panofsky  and  Dutton  and  Garratt  do  discuss  also 
the  so-called  Rossby  similarity  theory  summarized  above. 
The  derivations  assume  that  the  flow  in  the  surface  layer 
is  unidirectional  and  independent  of  /.  For  example,  the 
starting  point  for  McWilliams’  derivation  is  based  on  the 
idea  that  the  mean  vector  velocity  profile  u(z)  has  a  large 
shear  with  a  profile  shape  governed  by  the  boundary  stress 
(characterized  by  the  friction  velocity  it*)  and  the  near¬ 
boundary  turbulent  eddy  size  (effectively  the  height  z)  in 
the  following  way: 


du  it*  „ 
dz  KZ 


(9) 


where  s  is  a  unit  vector  in  the  direction  of  the  surface  shear 
stress  and  other  quantities  are  defined  above.  This  equation 
may  be  integrated  to  yield 


u(z)  =  —  log  (— )  s,  (10) 

K  \ZoJ 

where,  again,  z0  is  the  roughness  length.  Accordingly, 
the  wind  is  uni-directional  in  the  direction  of  the  surface 
shear  stress  and  increases  logarithmically  in  magnitude 
with  height.  Note  that,  if  expressed  as  wind  components  in 
any  locally  orthogonal  coordinate  system,  both  components 
must  increase  in  magnitude  with  height. 


4 If  the  mean  radial  flow  above  the  boundary  layer  were  not  zero,  the 
tangential  flow  would  evolve  with  time  on  account  of  the  material 
conservation  of  absolute  angular  momentum,  except  in  the  special  case 
of  a  vortex  in  which  the  flow  above  the  boundary  layer  is  along  absolute 
angular  momentum  surfaces. 
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Bulk  scheme  108-120  li 


(b)  0  20  40  60  80  100 


(d)  0  20  40  60  80  100 


Linear  model 


(f)  radius  (km) 


Figure  1.  Radius-height  cross-sections  of  azimuthally-averaged  radial  (thin/blue  contours)  and  tangential  (thick/red  contours)  wind  components  in  the 
lowest  2  km  averaged  at  15  minute  intervals  during  the  period  108-120  hours  for  the  different  boundary  layer  schemes,  (a)  bulk  scheme,  (b)  Blackadar 
scheme,  (c)  Burk-Thompson  scheme,  (d)  MRF  scheme,  (e)  Gayno-Seaman  scheme,  and  (f)  the  steady  linear  model.  Contour  interval  5  m  s"  1 . 


The  assumption  of  a  constant  (vector)  stress  throughout 
the  surface  layer  not  only  precludes  the  possibility  that 
the  wind  turns  with  height  in  that  layer,  it  precludes  also 
the  possibility  that  the  vertical  gradient  of  the  radial  and 
tangential  wind  components  can  have  different  signs,  or 
the  vertical  gradient  of  the  radial  component  can  be  zero. 
As  far  as  we  are  aware,  the  validity  of  this  assumption  in 
the  turbulent  boundary  layer  of  a  rapidly-rotating  vortex 
has  not  been  demonstrated  and  we  present  evidence  below 
from  various  numerical  calculations  as  well  as  observations 
suggesting  that  the  assumption  cannot  be  justified  in  the 
inner  core  of  a  hurricane. 

2.2.  Near-surface  wind  structure  in  mesoscale  models 

There  is  both  observational  and  theoretical  support  to 
suggest  that  the  vertical  gradients  of  the  radial  and 
tangential  wind  components  have  different  signs  near  the 
surface  in  the  inner  core  of  a  hurricane,  a  feature  that  is 
not  compatible  with  one  property  of  the  surface  layer  noted 
above.  The  tendency  to  produce  the  maximum  radial  inflow 
at  the  surface  is  evident  for  all  the  schemes  investigated  by 
Braun  and  Tao  (2000),  Nolan  et  al.  (2009a,b)  and  Smith 
and  Thomsen  (2010),  whether  or  not  the  formulation  of 
the  scheme  incorporated  a  log  layer.  As  an  illustration  of 
this  feature,  we  show  in  Figure  1  vertical  cross  sections 
of  the  azimuthally-averaged  radial  and  tangential  wind 
speed  components  in  the  idealized  hurricane  simulations 
described  by  Smith  and  Thomsen  op.  cit.  for  five  different 
boundary  layer  schemes.  The  cross  sections  encompass  the 
lowest  2  km  in  height  with  the  velocity  fields  being  averaged 
at  15  minute  intervals  during  the  mature  stage  of  vortex 
evolution  (the  period  108-120  hours).  The  boundary-layer 
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schemes  include  the  bulk  scheme,  the  Blackadar  scheme, 
the  Burk-Thompson  scheme,  the  MRF  scheme,  and  the 
Gayno-Seaman  scheme,  details  of  which  are  summarized  in 
Smith  and  Thomsen  with  references.  Kepert  (2012)  gives  an 
erudite  summary  of  the  essential  features  of  these  different 
schemes. 

For  comparison.  Figure  If  shows  an  example  of  cross 
sections  obtained  by  solving  the  quasi-linear  boundary- 
layer  model  with  a  prescribed  tangential  wind  profile5  just 
below  the  top  of  the  layer.  In  this  example,  the  maximum 
tangential  wind  speed  at  large  height  is  taken  to  be  60  m 
s_1  and  the  eddy  diffusivity  is  taken  to  be  a  constant,  equal 
to  100  m2  s_1.  The  quasi-linear  model  is  locally  analogous 
to  the  classical  Ekman  layer  model6  (Eliassen  and  Lystadt 
1977,  Kepert  2001).  Although  it  has  been  shown  that  the 
quasi-linear  approximation  becomes  invalid  in  the  inner 
core  of  a  hurricane  (Vogl  and  Smith  2009),  this  model  shows 
also  the  tendency  to  produce  the  maximum  radial  wind 
component  at  the  surface. 

As  an  aid  to  comparing  the  schemes  in  Figure  1, 
we  show  in  Figure  2  the  vertical  profiles  of  the  radial 


5Profile  3  in  Figure  1  of  Smith  (2003). 

f>The  quasi-linear  model  for  the  boundary  layer  is  obtained  by  neglecting 
the  nonlinear  acceleration  terms  for  the  agradient  wind  in  the  horizontal 
momentum  equations  (7)  and  (8)  and  the  centrifugal  and  Coriolis  terms 
are  linearized  about  the  gradient  wind  of  the  bulk  vortex  at  the  top 
of  the  boundary  layer  (Vogl  and  Smith  2009).  The  radial  diffusion  of 
momentum,  is  neglected  also.  These  equations  have  the  form  —C,av'  = 
(9 / dz)(Kdu/ dz)  and  —  £u  =  ( d / dz)( K dv'  /dz),  where  u  and  v'  are 
the  radial  and  tangential  components  of  the  ageostrophic  wind.  £a  is  the 
absolute  vorticity  of  the  gradient  wind,  £  =  2v/r  +  /  is  twice  the  absolute 
angular  velocity  of  the  gradient  wind,  K  is  the  vertical  eddy  diffusivity 
and  z  is  the  height.  Simple  closed-form  solutions  may  be  obtained  if  the 
diffusivity  is  assumed  to  be  constant  with  height. 
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and  tangential  wind  components  at  the  radius  ( rmax )  of 
the  maximum  azimuthally-averaged  tangential  wind  speed 
( Vmax )■  Included  also  are  the  corresponding  profiles  for 
the  quasi-linear  solution  shown  in  Figure  If.  Again  these 
profiles  highlight  the  fact  that  in  all  cases,  vmax  occurs  near 
the  top  of  the  inflow  layer.  It  is  particularly  noteworthy 
that  for  all  schemes,  and  for  the  quasi-linear  solution,  the 
maximum  radial  wind  speed  occurs  at  or  very  close  to  the 
surface. 


30  40  50  60  70  80  -40  -20  0  20 

v  (m/sec)  u  (m/sec) 

Figure  2.  Vertical  profiles  of  azimuthally-averaged  radial  and  tangential 
wind  components  in  the  lowest  2  km  at  the  radius  of  maximum  tangential 
wind  speed  for  the  different  panels  in  Figure  1.  (b)  bulk  scheme,  (B) 
Blackadar  scheme,  (BT)  Burk-Thompson  scheme,  (M)  MRF  scheme,  (GS) 
Gayno-Seaman  scheme  and  (L)  steady  quasi-linear  model. 

It  is  significant  that  for  the  region  inside  the  radius  of 
maximum  tangential  winds,  Kepert’s  (2012)  solutions  have 
also  the  tendency  to  produce  the  maximum  radial  inflow 
at  the  surface.  This  is  a  feature  of  all  the  schemes  he 
investigated,  whether  or  not  a  log  layer  was  “imposed” 
by  the  choice  of  the  linear  variation  of  near-surface  eddy 
diffusivity  with  height  (see  e.g.  his  Figures  3  to  6).  As  noted 
above,  such  a  feature  is  inconsistent  with  the  existence  of 
a  log  layer.  Even  so,  it  is  pertinent  to  mention  that  the 
numerical  model  used  by  Kepert,  as  well  as  the  quasi- 
linear  boundary-layer  model  summarized  above,  have  an 
issue  in  that,  as  pointed  out  by  Smith  and  Montgomery 
(2010),  they  effectively  prescribe  the  tangential  wind  at  the 
top  of  the  boundary  layer  where  the  flow  is  upwards.  Such 
a  prescription  at  an  outflow  boundary  makes  the  physical 
problem  ill-posed  as  the  boundary  layer,  itself,  should  be 
allowed  to  determine  the  tangential  momentum  that  it  expels 
into  the  bulk  vortex  aloft  (see  also  Rotunno  and  Bryan  2012, 
pl7). 

As  noted  above,  a  property  of  the  surface  layer  defined  by 
the  solution  (6)  is  the  strict  unidirectional  nature  of  the  wind 
within  it.  However,  if  the  solution  (5)  continues  to  hold  for 
some  distance  above  the  matching  layer,  the  wind  profile 
may  remain  logarithmic  with  height  while  accommodating 
some  cyclonic  turning  with  height  (Blackadar  and  Tennekes 
1968). 

To  show  that  the  unidirectional  surface  layer  is  not  a 
feature  of  any  of  the  parameterization  schemes  in  Figure 
1,  irrespective  of  whether  they  are  able,  in  principle,  to 
represent  a  log-layer,  we  show  in  Figure  3  the  hodographs 
of  the  wind  profiles  in  Figure  2.  Except  in  the  MRF-scheme 
(which,  as  noted  by  Smith  and  Thomsen,  predicts  a  rather 
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diffuse  boundary  layer)  and  in  the  quasi-linear  solution, 
the  schemes  indicate  that  a  significant  turning  of  the  wind 
vector  with  height  occurs  in  the  lowest  few  hundred  metres, 
a  property  that  cannot  be  represented  by  the  traditional 
surface  layer.  Note  that  except  in  panel  (c)  of  Figure  3,  the 
radial  wind  component  remains  approximately  constant  or 
actually  decreases  in  magnitude  with  height  in  the  lowest 
50  m.  This  feature  is  not  a  property  of  the  traditional  log- 
layer,  where  both  components  must  increase  with  height. 
Note  also  that  the  mean  wind  in  this  layer  is  generally  not 
in  the  direction  of  the  surface  wind,  which  again  is  not  a 
property  of  the  traditional  surface  log-layer. 

3.  Observations  of  the  hurricane  boundary  layer 

The  data  used  by  Powell  et  al.  (2003)  to  justify  the  presence 
of  a  log-layer  have  a  great  deal  of  scatter  and  individual 
profiles  are  not  shown.  However,  other  studies  indicate 
that  many  individual  inner-core  wind  soundings  do  not 
exhibit  the  structure  of  a  traditional  log-layer.  In  fact,  Global 
Positioning  System  (GPS)  dropwindsonde  data  in  the  inner 
core  of  hurricanes  and  typhoons  frequently  show  radial 
wind  profiles  that  have  a  maximum  inflow  at  the  surface. 

Composite  dropwindsonde  soundings  in  the  eyewall 
of  individual  storms  can  provide  a  useful  perspective 
on  the  vertical  structure  of  the  boundary  layer  winds 
there.  Two  such  examples  are  illustrated  here  using  data 
from  Typhoon  langmi  (Figure  4)  presented  by  Sanger 
et  al.  (2012),  and  Hurricane  Isabel  (2003)  (Figure  5) 
presented  by  Montgomery  et  al.  (2006)  and  Bell  and 
Montgomery  (2008).  In  both  cases,  the  maximum  storm- 
relative  tangential  wind  speed  (vmax)  occurs  within  the 
layer  of  relatively  strong  inflow  ( u ).  Without  exception,  the 
tangential  wind  component  ( v )  is  a  minimum  at  the  surface. 
While  the  magnitude  of  v  in  langmi  increases  with  height 
near  the  surface,  that  of  u  decreases  with  height,  except  in  a 
very  shallow  layer  (below  50  m)  in  the  supertyphoon  stage. 
In  Isabel,  a  negative  vertical  gradient  of  radial  velocity  is 
evident  throughout  the  boundary  layer  for  two  out  of  three 
days,  except  in  a  very  shallow  layer  below  50  m.  On  14 
September,  the  maximum  inflow  resides  at  the  surface. 

Figure  6  shows  the  hodographs  of  the  eyewall  wind 
composites  for  Jangmi  presented  in  Figure  4.  As  in 
the  numerical  calculations  shown  in  Figure  3,  much  of 
the  turning  of  the  wind  occurs  within  the  lowest  few 
hundred  metres.  The  turning  is  particularly  marked  in  the 
supertyphoon  stage  of  Jangmi,  for  which  the  hodograph  is 
quite  similar  to  that  in  the  bulk  scheme  (compare  Figure 
3a  with  Figure  6c).  This  finding  would  indicate  that  the 
bulk  scheme  is  not  necessarily  as  poor  as  Kepert  ’s  (2012) 
critique  of  it  might  suggest. 

The  foregoing  properties  are  reasons  to  question  the 
existence  of  a  surface-based  log-layer  in  which  the  wind  and 
shear  stress  vector  are  unidirectional. 

It  may  be  argued  that  Typhoon  Jangmi  and  Hurricane 
Isabel  are  only  two  storms  and  that  the  structures 
shown  in  Figures  4-6  may  not  be  general.  For  this 
reason  we  show  in  Figure  7  composite  plots  of  storm- 
relative  radial  and  tangential  velocity  for  the  eyewall 
region  of  thirteen  Atlantic  hurricanes.  These  eyewall 
composite  profiles  were  constructed  from  data  used  to 
characterize  the  mean  boundary  layer  structure  of  the 
near-core  vortex  region  contained  within  a  radius  of 
about  four  times  the  radius  of  maximum  tangential  winds 
(Zhang  et  al.  2011a).  The  eyewall  composites  consist  of 
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Bulk  scheme  Blackadar  scheme  Burk— Thompson  scheme 


Figure  3.  Wind  hodographs  in  the  lowest  2  km  corresponding  with  the  vertical  profiles  in  Figure  2.  (a)  bulk  scheme,  (b)  Blackadar  scheme,  (c)  Burk- 
Thompson  scheme,  (d)  MRF  scheme,  (e)  Gayno-Seaman  scheme  and  (f)  Quasi-linear  model.  The  tick  marks  on  the  curves  indicate  height  intervals  every 
50  m  stalling  at  the  surface  and  ending  at  400  m.  The  blue  and  green  lines  indicate  the  wind  vectors  at  the  surface  and  at  a  height  of  400  m,  respectively. 


v  (m/sec)  u  (m/sec) 

Figure  4.  Vertical  profiles  of  the  storm-relative  tangential  (v)  and  radial 
(u)  wind  components  in  the  eyewall  composites  in  Typhoon  Jangmi 
(2008).  Indices  1,2,3  on  the  curves  denote  the  tropical  storm,  typhoon  and 
supertyphoon  stages,  respectively.  Data  courtesy  of  N.  T.  Sanger. 
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Figure  5.  Vertical  profiles  of  the  storm-relative  tangential  (v)  and  radial 
(u)  wind  components  in  the  eyewall  composites  in  Hurricane  Isabel  (2003) 
on  three  consecutive  days  of  observations  (12-14  September).  Numbers  on 
curves  denote  the  date.  Data  courtesy  of  M.  A.  Bell. 


several  hundred  dropwindsondes.  In  these  composites,  the  height  near  the  surface.  The  radial  velocity  component 
tangential  velocity  component  increases  in  magnitude  with  increases  slightly  in  magnitude  with  height  within  the  first 
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50  metres,  and  subsequently  decreases  rapidly.  The  increase 
of  the  two  wind  components  in  the  lowest  50  metres  would 
not  rule  out  the  existence  of  a  log  layer  there.  To  examine 
this  possibility,  we  plot  in  Figure  8  the  total  wind  speed 
from  these  components  as  a  function  of  height  in  the  lowest 
kilometre  and  also  the  wind  hodograph  to  a  height  of 
two  kilometres.  It  is  seen  that  while  the  wind  profile  is 
approximately  logarithmic  in  the  layer  between  100  and  400 
m,  this  logarithmic  profile  does  not  extend  all  the  way  to  the 
surface.  Moreover,  as  in  the  hodographs  shown  in  Figures  3 
and  6,  the  wind  vector  turns  through  an  appreciable  angle 
within  this  layer,  ruling  out  that  the  layer  behaves  as  a 
constant  stress  layer.  While  it  might  be  argued  that  the 
logarithmic  behaviour  would  be  consistent  with  solution 
(5),  the  decrease  in  the  radial  component  of  flow  above  50 
m  is  evidence  that  horizontal  inhomogeneities  associated 
with  the  decrease  in  the  agradient  wind  with  height  are 
important. 

In  summary,  the  foregoing  observations  indicate  a 
significant  turning  of  the  wind  vector  with  height  in  the 
lowest  few  hundred  metres  of  the  inner-core  boundary 
layer,  generally  accompanied  by  a  decrease  in  the  radial 
wind  component  with  height.  These  features,  which  support 
the  modeling  results  discussed  in  section  2,  cannot  be 
represented  by  the  traditional  surface  log-layer.  Indeed,  in 
the  eyewall  composite  for  many  hurricanes,  the  vertical 
profile  of  storm-relative  wind  speed  does  not  strictly  follow 
a  logarithmic  profile  throughout  the  surface  layer. 

4.  Ramifications  for  hurricane  modelling 

Kepert  (2012)  criticizes  the  bulk  and  high  res  schemes  on 
one  ground  that  they  do  not  produce  the  observed  near¬ 
surface  logarithmic  layer.  As  argued  above,  the  basis  for  a 
log  layer  in  the  inner  core  of  a  hurricane  is  not  compelling, 
either  on  theoretical  or  observational  grounds.  Even  so,  the 
question  remains:  how  important  are  the  structural  details 
of  the  shallow  surface  layer  on  the  prediction  of  vortex 
evolution,  provided  that  the  surface  stress  and  surface  heat 
fluxes  are  adequately  represented? 

From  an  elementary  perspective,  in  the  classical  Ekman 
solution  (including  the  version  where  the  surface  stress  is 
assumed  to  be  in  the  direction  of  the  surface  wind),  the 
radial  volume  transport  depends  only  on  the  surface  stress 
and  not  on  the  details  of  the  shallow  surface  layer.  This 
result  follows  directly  by  integrating  the  steady  linearized 
tangential  momentum  equation  with  respect  to  height  and 
assuming  that  the  tangential  wind  approaches  the  gradient 
wind  at  large  heights.  A  similar  result  is  true  for  the  quasi- 
linear  model  of  the  boundary  layer  discussed  above. 

In  the  classical  Ekman  solution,  the  vertical  velocity  at 
the  top  of  the  boundary  layer  is  simply  proportional  to  the 
radial  gradient  of  the  volume  transport  in  the  layer,  and 
hence  to  the  radial  gradient  of  the  surface  shear  stress.  This 
result  follows  immediately  by  integrating  the  continuity 
equation  with  respect  to  height  (Gill  1982,  section  9.4).  Like 
the  volume  transport,  the  vertical  velocity  at  large  height 
does  not  depend  on  the  details  of  the  surface  layer.  These 
same  remarks  apply  also  to  the  quasi-linear  vortex  boundary 
layer  model  discussed  in  section  2. 

While  a  scale  analysis  shows  that  neither  the  Ekman 
model  or  the  quasi-linear  model  are  valid  in  the  inner  core 
region  of  a  hurricane  and  that  the  non-linear  acceleration 
terms  in  the  boundary-layer  equations  are  important  in  this 
region  (Smith  1968,  Vogl  and  Smith  2009),  it  has  yet  to  be 
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demonstrated  that  the  details  of  the  shallow  surface  layer 
have  a  profound  effect  on  the  volume  of  air  converging  in 
the  boundary  layer  and  hence  on  the  vertical  motion  out 
of  the  boundary  layer,  provided  that  the  surface  stress  is 
represented  appropriately. 

Of  course,  the  magnitude  and  vertical  distribution  of  the 
radial  velocity  component  depends  in  part  on  the  assumed 
vertical  profile  of  diffusivity,  even  in  the  classical  Ekman 
solution.  In  particular,  the  bulk  magnitude  of  the  diffusivity 
together  with  the  Coriolis  parameter  determines  the  depth 
of  the  inflow  layer  (e.g.  Gill  1982)  and  hence  the  depth 
over  which  the  volume  flux  is  distributed.  This  dependence 
on  diffusivity  extends  therefore  to  the  radial  advection  of 
absolute  angular  momentum.  Clearly,  the  efficacy  of  the 
boundary-layer  spin-up  mechanism  articulated  by  Smith 
et  al.  (2009)  will  depend  quantitatively  on  the  particular 
parameterization  scheme  as  confirmed  by  the  calculations 
of  Smith  and  Thomsen  (2010)  summarized  in  Figures  1-3 
herein. 

As  noted  above,  the  calculations  of  Smith  and  Thomsen 
(2010)  showed  that  a  recently  articulated  boundary-layer 
spin-up  mechanism  for  the  hurricane  by  Smith  et  al.  (2009) 
transcends  the  presence  of  a  log  layer.  Independently,  Bryan 
(2012,  his  Fig.  16)  has  shown  that  the  incorporation  of  a 
reduced  vertical  mixing  length  near  the  surface  (using  the 
Blackadar  formulation  for  vertical  mixing  length)  yields 
an  essentially  similar  dependence  of  maximum  tangential 
winds  on  the  ratio  of  enthalpy  and  momentum  surface 
exchange  coefficients,  while  the  simulated  vmax  tended 
“to  be  slightly  lower  with  the  Blackadar  formulation  for 
vertical  mixing  length.”  Together,  these  results  suggest  that 
the  essence  of  tropical  cyclone  spin  up  and  the  dependence 
of  maximum  winds  on  the  ratio  of  enthalpy/momentum 
exchange  coefficients  is  captured  without  a  log  layer. 

It  is  evident  from  the  foregoing  discussion  that 
uncertainties  in  the  optimum  scheme  for  use  in  operational 
hurricane  models  remain.  In  general,  such  models  do  not 
have  the  vertical  grid  resolution  to  resolve  the  log  layer, 
which  is  used  merely  to  extrapolate  the  wind  from  the 
lowest  model  grid  level  to  a  standard  height  (normally  10  m) 
where  the  normal  aerodynamic  drag  formulation  (in  terms 
of  Cd)  is  applicable.  In  view  of  the  large  uncertainties  in  the 
formulation  of  the  eddy  diffusivity  above  the  lowest  model 
level,  the  possible  sensitivity  to  the  precise  formulation 
of  the  surface  layer  may  be  over-exaggerated.  Indeed,  a 
major  issue  confronting  hurricane  modellers  is  the  lack  of 
observational  data  on  the  radial  and  vertical  structure  of 
eddy  diffusivity  in  the  strong  wind  region  of  hurricanes  to 
guide  the  formulation  of  this  quantity  in  models.  One  step  in 
this  direction  has  been  taken  in  two  recent  papers  by  Zhang 
et  al.  (201  lb)  and  Zhang  and  Montgomery  (2012). 

Despite  Kepert’s  (2012)  critique  of  the  bulk  and 
Blackadar  schemes,  the  results  of  section  2  show  that  the 
predictions  using  the  Blackadar  scheme  do  not  deviate 
significantly  from  the  other  schemes  and  even  those  of 
the  bulk  scheme  are  not  totally  unrealistic  compared  with 
some  of  the  observations  shown  above.  In  particular,  we 
have  shown  that  the  wind  speed  profiles  and  hodographs 
in  the  vicinity  of  the  eyewall  region  using  the  bulk  and 
Blackadar  schemes  are  not  inconsistent  in  magnitude  with 
those  in  major  hurricanes  such  as  Isabel  and  Jangmi.  While 
we  do  not  wish  to  defend  the  use  of  simple  boundary- 
layer  schemes  for  their  accuracy  in  operational  prediction 
models,  we  do  believe  that  they  have  an  important  role 
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Figure  6.  Wind  hodographs  in  the  lowest  2  km  corresponding  with  the  vertical  wind  profiles  of  Jangmi  shown  in  Figure  4.  The  tick  marks  on  the  curves 
indicate  height  intervals  every  50  m  starting  at  the  surface  and  ending  at  400  m.  The  blue  and  green  lines  indicate  the  wind  vectors  at  the  surface  and  at 
a  height  of  400  m,  respectively. 


in  generating  understanding  of  tropical-cyclone  spin  up 
and  maximum  potential  intensity.  We  believe  that  Kepert’s 
erudite  comparison  of  the  different  schemes  is  an  important 
step  in  attempts  to  determine  an  optimum  scheme  for  use  in 
operational  prediction  models.  Nevertheless,  for  the  reasons 
articulated  herein,  we  do  not  subscribe  to  his  assertion  that 
the  absence  of  a  log  layer  should  be  a  criterion  for  rejecting 
a  scheme. 

Another  issue  raised  by  our  results  is  the  validity  of 
assuming  a  constant  stress  layer  with  a  logarithmic  wind 
speed  profile  for  estimating  the  drag  coefficient  at  major 
hurricane  wind  speeds  (e.g.  Powell  2003,  Holthuijsen  et  al. 
2012).  The  basis  of  this  assumption  is  that  the  flow  can  be 
treated  as  horizontally  homogeneous,  which  we  have  shown 
here  to  be  untenable  on  both  observational  and  theoretical 
grounds.  In  particular,  the  nonlinear  inertial  effects  are 
shown  to  be  important  near  the  surface,  where  the  effective 
radial  pressure  gradient  force  is  largest.  A  consequence  is 
that  the  radial  flow  tends  to  be  a  maximum  at  or  near 
the  surface  as  seen  in  the  observations  and  the  vertical 
gradients  of  the  magnitude  of  the  radial  and  tangential  wind 
components  tend  to  have  opposite  signs.  A  method  that 
avoids  these  assumptions  in  estimating  the  drag  coefficient 
at  high  wind  speeds  is  discussed  by  Bell  et  al.  (2012), 
although  this  method  has  its  own  limitations  as  well.  This 
method  is  based  on  a  control  volume  analysis  of  absolute 
angular  momentum  and  total  energy  around  the  eyewall 
region  in  the  lower  troposphere. 

5.  Conclusions 

We  have  questioned  the  validity  of  the  traditional 
logarithmic  surface  layer  in  the  inner  core  of  hurricanes. 
Definitive  observational  evidence  for  its  existence  in 
previous  studies  is  tenuous  and  is  based  on  data  that  have 
a  significant  amount  of  scatter.  Indeed,  many  individual 
eyewall  soundings  and  a  composite  comprising  thirteen 
Atlantic  hurricanes  do  not  support  its  existence.  There 
are  theoretical  reasons  why  the  logarithmic  layer  may  be 
violated  in  the  inner  core  of  hurricanes:  this  is  because  the 
inward-directed  effective  pressure  gradient  force  is  largest 
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at  the  surface,  where  the  tangential  wind  is  reduced  the  most 
from  its  gradient  value  aloft.  The  existence  of  this  cross¬ 
stream  force  raises  the  possibility  that  the  largest  inflow 
occurs  at,  or  very  close  to,  the  surface,  which  would  imply 
that  the  horizontal  shear-stress  vector  is  not  unidirectional 
near  the  surface.  We  have  presented  both  numerical  model 
results  and  observational  analyses  in  support  of  these 
ideas.  We  noted  that  deviations  from  a  logarithmic  layer 
in  the  inner  core  of  hurricanes  described  herein  must 
affect  the  ability  to  infer  the  surface  drag  coefficient  from 
dropwindsonde  wind  profiles  using  methods  that  assume  a 
logarithmic  layer  from  the  outset.  Finally,  we  drew  attention 
to  a  study  examining  a  range  of  boundary-layer  schemes 
demonstrating  that  a  recently  articulated  boundary-layer 
spin-up  mechanism  transcends  the  presence  of  a  logarithmic 
layer. 


30  40  50  60  -30  -20  -iO  0  10 

v  (m/sec)  u  (m/sec) 

Figure  7.  Vertical  profiles  of  the  tangential  (v)  and  radial  (u)  wind 
components  in  the  eyewall  composites  of  many  hurricanes.  Data  courtesy 
J.  A.  Zhang. 
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Figure  8.  (a)  Vertical  profiles  of  the  tangential  (v)  and  radial  (u)  wind 
components  shown  in  Figure  7,  but  plotted  on  a  logarithmic  vertical  scale, 
(b)  Hodograph  of  the  wind  components  plotted  in  Figure  7.  Horizontal  tick 
marks  indicate  heights  from  10  m  to  400  m  at  10  m  intervals. 
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